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Abstract. We investigate how properties of the corona and solar wind in open coronal holes 
depend on properties of magnetic fields and their footpoint motions at the surface. We per- 
form one-dimensional magnetohydrodynamical (MHD) simulations for the heating and the ac- 
celeration in coronal holes by low-frequency Alfven waves from the photosphere to 0.3 or 0.1 AU. 
We impose low-frequency ( ~ 0.05Hz) transverse fluctuations of the field lines at the pho- 
tosphere with various amplitude, spectrum, and polarization in the open flux tubes with dif- 
ferent photospheric field strength, B r ,o, and super-radial expansion of the cross section, / max - 
We find that transonic solar winds are universal consequences. The atmosphere is also stably 
heated up to ~ 10 6 K by the dissipation of the Alfven waves through compressive-wave 
generation and wave reflection in the cases of the sufficient wave input with photospheric am- 
plitude, {dv±fi) ~ 0.7km s _1 . The density, and accordingly the mass flux, of solar winds 
show a quite sensitive dependence on (dv±fl) because of an unstable aspect of the heating 
by the nonlinear Alfven waves. A case with (dv±,o) = 0.4km s _1 gives ~50 times smaller 
mass flux than the fiducial case for the fast wind with (dv± t o) = 0.7km s _1 ; solar wind vir- 
tually disappears only if (dv±,o) becomes ~l/2. We also find that the solar wind speed has 
a positive correlation with B r ,o//max> which is consistent with recent observations by Kojima 
et al. Based on these findings, we show that both fast and slow solar winds can be explained 
by the single process, the dissipation of the low-frequency Alfven waves, with different sets 
of (duj_,o) and B r fi/f max . Our simulations naturally explain the observed (i) anti-correlation 
of the solar wind speed and the coronal temperature and (ii) larger amplitude of Alfvenic fluc- 
tuations in the fast wind. In Appendix, we also explain our implementation of the outgoing 
boundary condition of the MHD waves with some numerical tests. 



1. Introduction 

The corona can be roughly divided into two sectors with respect 
to magnetic field configurations on the solar surface. One is closed 
field regions in which both footpoints of each field line are rooted 
at the photosphere. The other is open field regions in which one 
footpoint is anchored at the photosphere and the other is open into 
the interplanetary space. The open field regions roughly coincide 
with areas called coronal holes which are dark in X-rays. From the 
coronal holes, the hot plasma streams out as solar winds. 

The Alfven wave, generated by the granulations or other sur- 
face activities, is a promising candidate operating in the heating 
and acceleration of solar winds from coronal holes. It can travel 
a long distance so that the dissipation plays a role in the heating 
of the solar wind plasma as well as the lower coronal plasma, in 
contrast to other processes, such as magnetic reconnection events 
and compressive waves, the heating of which probably concentrates 
at lower altitude. High-frequency (~ 10 4 Hz) ioncyclotron waves 
were recently suggested for the inferred preferential heating of the 
perpendicular temperature of the heavy ions [Axford & McKenzie 
, 1997; Kohl et al., 1998], whereas there is still a possibility that the 
observed temperature anisotropy might be due to a line-of-sight ef- 
fect [Raouafi & Solanki , 2004]. However, the protons which com- 
pose the main part of the plasma are supposed to be mainly heated 
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by low-frequency ( ~ 0.1Hz) components in the magnetohydro- 
dynamical (MHD) regime. This is first because the low-frequency 
wave is expected to have more power, and second because the res- 
onance frequency of the proton is higher than those of heavier ions 
so that the energy of the ioncyclotron wave is in advance absorbed 
by the heavy ions [Cranmer , 2000]. 

The treatment of the wave dissipation is quite complicated even 
for MHD processes. The main reason is that the nonlinearity 
becomes essential because amplitudes of upgoing waves are in- 
evitably amplified in the stratified atmosphere with decreasing den- 
sity. Therefore, numerical simulation is a powerful tool; there 
are several works on the wave propagation and dissipation in one- 
dimensional (ID) [Lau & Sireger , 1996; Kudoh & Shibata , 1999; 
Orta et al., 2003] and in two-dimensional [Ofman & Davila , 
1997, 1998; Ofman , 2004] simulations. However, large density 
contrast amounting more than 15 orders of magnitude from the 
photosphere to the outer heliosphere has prevented one so far from 
carrying out numerical simulation in a broad region from the pho- 
tosphere to the interplanetary space even in ID modeling. 

In Suzuki & Inutsuka [2005] (paper I hereafter) we successfully 
carried out ID MHD simulation with radiative cooling and thermal 
conduction from the photosphere to 0.3 AU. We showed that the 
coronal heating and the fast solar wind acceleration in the coronal 
hole are natural consequences of the transverse footpoint fluctua- 
tions of the magnetic field lines. Low-frequency ( ~ 0.05Hz) 
Alfven waves are generated by the photospheric fluctuations [Ul- 
rich , 1996; Cranmer & van Ballegooijen , 2005]. The sufficient 
amount of the energy transmits into the corona in spite of attenu- 
ation in the chromosphere and the transition region (TR), so that 
they effectively dissipate to heat and accelerate the coronal plasma. 
Although it is based on the one-fluid MHD approximation, it, for 
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the first time, self-consistently treats the plasma heating and the 
propagation of the Alfven waves from the photosphere to the inter- 
planetary space. However, we studied only one case for the fast so- 
lar wind in paper I. It is important to examine how properties of the 
corona and the solar wind are affected by adopting different types 
of the wave injections and/or different magnetic field geometry of 
the flux tubes. The purpose of this paper is to study parametrically 
the heating by the low-frequency Alfven waves. 

In this paper we adopt the one-fluid MHD approximation in ID 
flux tubes, which is the same as in paper I. This assumption is not 
strictly correct in the solar wind plasma, which will be discussed 
later (§6.2). However, we think that the MHD approximation is ap- 
propriate for studies of the average properties of the plasma in the 
corona and the solar wind because the plasma particles are suffi- 
ciently randomized by the fluctuating magnetic fields. 

2. Simulation Set-up 

We consider ID open flux tubes which are super-radially open, 
measured by heliocentric distance, r. The simulation regions are 
from the photosphere (r = IRq) with density, p = lCF 7 g cm -3 , 
to 657? (0.3AU) or ~ 20R Q (~0.1AU), where R Q is the solar 
radius. Radial field strength, B r , is given by conservation of mag- 
netic flux as 



B r r f(r) = const., 



(1) 



where f(r) is a super-radial expansion factor. We adopt the same 
function as in Kopp & Holzer [1976] for f(r) but consider two-step 
expansions: 



fir) = 



/i,max exp( 2 



+ /i./Wexp(^) + / 2 



exp(^) + 1 



exp(^) + 1 



where /i = 1 - (/i, max - 1) exp( R@ ai Rl ) and f 2 = 1 - 

(/2,max — 1) exp( R@ a — ). The flux tube initially expands by a 
factor of /i, m ax around Ri = 1.01.R© corresponding to the 'fun- 
nel' structure [Tu et al., 2005], and followed by /2, ma x times ex- 
pansion around R 2 = 1.2R Q due to the large scale dipole mag- 
netic fields (Figure 1). We define the total super-radial expansion 

aS /max — fl ,max/2,max • 

We prescribe the transverse fluctuations of the field line by the 
granulations at the photosphere, which excite Alfven waves. We 
study cases of various spectra, pol arization s, and root mean squared 
(rms) amplitudes, (dt>x,o)(= y/Jdv± ))(km s _1 ). (dv±,o) and 



power spectrum, P(u)(erg g x Hz x ), have a relation of 



P(v)dv. 



(2) 



The parameters adopted in different simulation Runs are tabulated 
in Table 1. We consider three types of spectrum : (i) P(y) oc v~ x 
(ii) P{y) oc v° (white noise) (iii) P(u) oc 5(1/ v - 180s) (si- 
nusoidal waves with a period of 180 seconds) . For the first and 
second cases we set v\ OVJ cut = 6 x 10~ 4 Hz (a period of ~ 28min) 
and v up cut = 0.05Hz (20s) for Run I or ^ up cut = 0.025Hz 
(40s) for the others. For the linearly polarized perturbations we 
only take into account one transverse component besides the ra- 
dial component, while we consider the two transverse components 
for the circularly polarized case. Values of (dvx.o) are chosen to 
be compatible with the observed photospheric velocity amplitude 
~ 1km s _1 [Holweger et al., 1978]. At the outer boundaries, non- 
reflecting condition is imposed for all the MHD waves (see Ap- 
pendix), which enables us to carry out simulations for a long time 
until quasi-steady state solutions are obtained without unphysical 
wave reflection. 

We dynamically treat the propagation and dissipation of the 
waves and the heating and acceleration of the plasma by solving 
ideal MHD equations with the relevant physical processes (paper 
I): 



dp . P d 2 
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(r 2 fF c ) + q R = 0, 







[r^/f(v±B r -v r B ± )], 



(6) 



(7) 



dt r*/] dr 

where p, v, p, B are density, velocity, pressure, and magnetic field 
strength, respectively, and subscripts r and _L denote radial and 
tangential components; and denote Lagrangian and Eulerian 



Table 1. Model Parameters 



Run 


I 


II 


III 


IV 


V 


Outer Bound. 


65R Q 


2OR 


20Rq 


2OR 


2OR 


Spectrum 


u- 1 


u- 1 


u° 


Sii/- 1 - 180s) 


S(u- X -180s) 


Polarization 


Linear 


Linear 


Linear 


Linear 


Circular 


(dv± o)(km/s) 


0.7 


0.7 


0.7 


0.7 


0.7 


Br,o(G) 


161 


161 


161 


161 


161 


/l.max 


30 


30 


30 


30 


30 


J2,max 


2.5 


2.5 


2.5 


30 


2.5 


Description 


Fast Wind 
(paper I) 


Coarse Ver. 
of Run I 









VI 


VII 


VIII 


IX 


X 


2Oi? 
v- 1 


2O,R 
u- 1 


22Rq 
u- 1 


2OR 


65R Q 
u- 1 


Linear 


Linear 


Linear 


Linear 


Linear 


1.4 


0.4 


0.7 


0.3 


1.0 


161 


161 


322 


161 


322 


30 


30 


45 


30 


45 


2.5 


2.5 


10 


2.5 


10 








No Steady 


Slow Wind 








State 
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derivatives, respectively; e = — — j- - is specific energy and we as- 
sume the equation of state for icfeal gas with a ratio of specific heat, 
7 = 5/3; G and M© are the gravitational constant and the solar 
mass; F c {= kqT 5 ^ 2 ^) is thermal conductive flux by Coulomb 
collisions, where ko = 1CT 6 in c.g.s unit [Braginskii , 1965]; </r is 
radiative cooling described below. We use optically thin radiative 
loss[Landini & Monsignori-Fossi, 1990] in the corona and upper 
TR where temperature, T > 4 x 10 4 K. In the chromosphere and 
low TR, we adopt empirical radiative cooling based on the obser- 
vations [Anderson & Athay , 1989; Moriyasu et al., 2004] to take 
into account the optically thick effect. In the low chromosphere, 
the temperature sometimes drops to ~ 3000K because we do not 
consider other heating sources, such as sound waves, which must 
be important there besides Alfven waves [Carlsson & Stein , 1992; 
Bogdan et al, 2003]. Such unrealistically low temperature inter- 
rupts the propagation of the Alfven waves in the low chromosphere 
through the change of the density structure. To give realistic esti- 
mates of the transmission of the Alfven waves there, we switch off 
the radiative cooling, if the temperature becomes < 5000K only 
when p > 10 _1 g cm" . 

We adopt the second-order MHD-Godunov-MOCCT scheme 
(Sano & Inutsuka 2006) to update the physical quantities. We solve 
eqs. (3)-(6) on fixed Eulerian mesh by remapping the physical vari- 
ables updated in Lagrangian coordinate onto the original grids at 
each time step. The induction equation (7) is solved on the Eule- 
rian grids. We use implicit time steps for the thermal conduction 
and radiative cooling in energy equation (6) because the Courant 
conditions are severe for these processes, while explicit time steps 
are adopted for the other terms. Each cell boundary is treated as 
discontinuity, and for the time evolution we solve nonlinear Rie- 
mann shock tube problems with the magnetic pressure term by us- 
ing the Rankin-Hugoniot relations. Therefore, entropy generation, 
namely heating, is automatically calculated from the shock jump 
condition. A great advantage of our code is that no artificial viscos- 
ity is required even for strong MHD shocks; numerical diffusion is 
suppressed to the minimum level for adopted numerical resolution. 

We initially set static atmosphere with a temperature T = 10 4 K 
to see whether the atmosphere is heated up to coronal temperature 
and accelerated to accomplish the transonic flow. At t = we start 
the inject of the transverse fluctuations from the photosphere and 
continue the simulations until the quasi-steady states are achieved. 

3. Heating and Acceleration by Wave 
Dissipation 

Before showing results of the various coronal holes, we explain 
how the coronal heating and the solar wind acceleration were ac- 
complished in the fiducial case (Run I) which was studied in pa- 
per I for the fast solar wind. Figure 2 plots the initial condition 
(dashed lines) and the results after the quasi-steady state condition 
is achieved at t = 2573 minutes (solid lines) in Run I, compared 
with recent observations of fast solar winds. From top to bottom, 
tv(km s _1 ), T(K), mass density, p(g cm -3 ), and rms transverse 
amplitude, (cfox)(km s _1 ) are plotted. As for the density, we com- 
pare our result with observed electron density, N e , in the corona. 
When deriving N e from p in the corona, we assume H and He 
are fully ionized, and 7V e (cm~ 3 ) = 6 x 10 23 p(g cm -3 ). These 
variables are averaged for 3 minutes to incorporate observational 
exposure time. 

Figure 2 shows that the initially cool and static atmosphere 
is effectively heated and accelerated by the dissipation of the 
Alfven waves. The sharp TR which divides the cool chromosphere 
with T ~ 10 4 K and the hot corona with T ~ 10 6 K is formed 
owing to a thermally unstable region around T ~ 10 5 K in the ra- 
diative cooling function [Landini & Monsignori-Fossi, 1990]. The 
hot corona streams out as the transonic solar wind. The simulation 
naturally explains the observed trend quite well, (see paper I for 
more detailed discussions.) 

Figure 3 shows transfer of the energy (upper) and momentum 
(lower) flux of various components described below. The energy 
equation in an Eulerian form becomes 



dE „ 



where E and F are the total energy density and flux. — Li oss is loss 
of the energy. For the analyses of quasi-steady state behaviors, we 
can reasonably assume that the time-average of the first term equals 
to 0. Then, in our case the explicit form becomes 



2 f dr 



2 j. 

r f 



pV r 



2 i 2 
V r +V, 



P 

+ e + 1 -- 



GMr, 



+ Bi Vr _B I B^ + Fc 

4-K 47T 



+ (qn) = 0, 



(8) 



where (■ ■ •) denotes the time-average. Using the relation, v± = 
—B^/^/Airp, for the outgoing Alfven wave, we can extract the 
energy flux of the Alfven wave from eq.(8) as follows [lacques , 
1977]: 



(Fa) = (pV 2 ± (vA + V r ) +PAV r ) = 



B r B±v± 
4tt 



B' 



A 2 8n 



Vr 



B\ 
8n 



(9) 



where we define Alfven wave pressure as pa = Bj_/8n. In the up- 
per panel of Figure 3 we plot the wave energy flux, {Fa) ^f(r\ ' 
normalized by cross section of the flux tube at r c = 1.02J? o 
with kinetic energy flux, pv r \ ^ , enthalpy flux, pv r (e + 



P\ r f(r) 
P> rif(r c ) 



F + Ll oss = 0, 



2 r'if(r c ) ' 

, and integrated radiative loss, /T q r r .> // r \ dr. Note 
that the real energy flux is larger (smaller) in r < r c (r > r c ). 
While the most of the initial wave energy escapes by radiation 
loss in the chromosphere and the low corona, the remained energy, 
which is ~ 10% of the input, is transfered to the kinetic energy of 
the solar wind in the outer region. The thermal energy is almost 
constant in the corona, which indicates that the ~1MK plasma is 
maintained by the energy balance among wave heating, radiative 
cooling, adiabatic cooling (solar wind acceleration), and redistri- 
bution by thermal conduction. 

The lower panel of Figure 3 plots measure of pressure of 
Alfven waves, {pAV r ), and thermal gas, (pv r ). This shows that 
the Alfven wave pressure dominates the gas pressure in the solar 
wind acceleration region (1.5J? ~ r ~ WRq); the fast so- 
lar wind is driven by the wave pressure rather than by the thermal 
pressure. 

In paper I, we claimed that the dissipation of the low-frequency 
Alfven waves in the corona is mainly by the generation of slow 
waves and reflection of the Alfven waves due to the density fluc- 
tuation of the slow modes. This can be seen in r — t diagrams. 
Figure 4 presents contours of amplitude of v r , p, v±, and B±/B r 
in R Q < r < 15R® from t = 2570 min. to 2600 min. Dark 
(light) shaded regions denote positive (negative) amplitude. Above 
the panels, we indicate the directions of the local 5 characteristics, 
two Alfven , two slow, and one entropy waves at the respective 
positions. Note that the fast MHD and Alfven modes degenerate 
in our case (wave vector and underlying magnetic field are in the 
same direction), so we simply call the innermost and outermost 
waves Alfven modes. In our simple ID geometry, v r and p trace 
the slow modes which have longitudinal wave components, while 
v± and B± trace the Alfven modes which are transverse (see Cho 
& Lazarian [2002;2003] for more general cases.). 

One can clearly see the Alfven waves in v± and B±/B r dia- 
grams, which have the same slopes with the Alfven characteristics 
shown above. One can also find the incoming modes propagating 
from lower-right to upper-left as well as the outgoing modes gener- 
ated from the surface 1 . These incoming waves are generated by the 
reflection at the 'density mirrors' of the slow modes (see paper I). 
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At intersection points of the outgoing and incoming characteristics 
the non-linear wave-wave interactions take place, which play a role 
in the wave dissipation. 

The slow modes are seen in v r and p diagrams. Although it 
might be difficult to distinguish, the most of the patterns are due 
to the outgoing slow modes 2 which are generated from the pertur- 
bations of the Alfven wave pressure, B\/%n [Kudoh & Shibata , 
1999; Tsurutani et al., 2002]. These slow waves steepen eventually 
and lead to the shock dissipation. 

The processes discussed here are the combination of the di- 
rect mode conversion to the compressive waves and the parametric 
decay instability due to three-wave (outgoing Alfven , incoming 
Alfven , and outgoing slow waves) interactions (Goldstein 1978; 
Terasawa et al. 1986; see also §4.1.2) of the Alfven waves. Al- 
though they are not generally efficient in the homogeneous back- 
ground since they are the nonlinear mechanisms, the density gra- 
dient of the background plasma totally changes the situation. Ow- 
ing to the gravity, the density rapidly decreases in the corona as 
r increases, which results in the amplification of the wave am- 
plitude so that the waves easily become nonlinear. Furthermore, 
the Alfven speed varies a lot due to the change of the density 
even within one wavelength of Alfven waves with periods of min- 
utes or longer. This leads to both variation of the wave pres- 
sure in one wavelength and partial reflection through the deforma- 
tion of the wave shape [Moore et al., 1991]. The dissipation is 
greatly enhanced by the density stratification, in comparison with 
the case of the homogeneous background. Thus, the low-frequency 
Alfven waves are effectively dissipated, which results in the heat- 
ing and acceleration of the coronal plasma. 

In summary, we have shown that the low-frequency Alfven waves 
are dissipated via the combination of the direct mode conversion 
and the decay instability which are considerably enhanced by the 
long wavelength (so called "non-WKB") effect in the stratified 
atmosphere. Although each process is not individually effective 
enough, the significant dissipation rate is realized by the combina- 
tion. 

4. Universality and Diversity 

We investigate how the properties of the corona and solar wind 
depend on the parameters of the input photospheric fluctuations and 
the flux tube geometry. We use smaller simulation boxes in Run 
II-IX than in Run I (Table 1) to save the computational time for 
the parameter studies. In most cases quasi-steady state conditions 
are achieved after t ~ 30hr. In this section we show results at 
t = 48.3hr unless explicitly stated. 

4.1. Various Input Waves 
4.1.1. Wave Spectrum 

We show dependences on the spectrum, P(u), of the photo- 
spheric fluctuations. We here compare the results of P(y) oc v~ x 
(6 x 10~ 4 < v < 2.5 x 10~ 2 Hz; Run II), //(white noise; Run 
III), and purely sinusoidal perturbation with a period of 180 s, 
P{v) oc — 180s) (Run IV). Figure 5 shows structure of 

the coronae and solar winds. The upper panel of Figure 6 plots an 
adiabatic constant, S c , of the outgoing Alfven wave derived from 
wave action [Jacques , 1977]: 



Figure 5 shows that the plasma is heated up to the coronal tem- 
perature and accelerated to 400-600km/s at 0.1 AU in all the Runs. 



Particularly, the results of P(y) 



and v are quite simi- 



2 f(r 



va r%f{r c ) ' 



where 



(10) 



(11) 



is amplitude of the outgoing Alfven wave (Elsasser variables). 
The lower panel shows nonlinearity, (Sva,+)/va of the outgoing 
Alfven waves. 



lar. The attenuation of the outgoing Alfven waves is due to the 
two kinds of the processes, the shocks and the reflection in our 
simulations. The shock dissipation per unit distance is larger for 
waves with larger v because the number of shocks is proportional 
to wavenumber. On the other hand, Alfven waves with smaller v 
suffer the reflection more since the wavelengths are relatively larger 
compared to the variation scale of the Alfven speed so that the wave 
shapes are deformed more. These two ingredients with respect to 
the ^-dependence cancel in these two cases, and the variations of 
the energy and the amplitude become quite similar (Figure 6). The 
structures of the corona and the solar wind are also similar since 
they are consequences of the wave dissipation. 

The results of the purely sinusoidal fluctuation (Run IV) are 
slightly different from the other two cases. The TR is not so sharp 
as those in the other two Runs. The temperature firstly rises at 
r — R Q — 0.005i? Q , but drops slightly and again rises around 
r — Rq — O.O2i? . Actually, the TR moves back and forth 
more compared to the other cases. This is because the shock heat- 
ing takes place rather intermittently only at the wave crests of the 
monochromatic wave, while in the other cases the heating is dis- 
tributed more uniformly by the contributions from waves with var- 
ious v's. As a result, the chromospheric evaporation occurs inter- 
mittently so that the TR becomes "dynamical". The sharpness of 
the TR is important in terms of the wave reflection. In the sinu- 
soidal case (Run IV) the Alfven waves becomes more transparent 
at the TR since the variation of the Alfven speed is more gradual. 
Therefore, more wave energy transmits into the corona, avoiding 
the reflection at the TR (Figure 6). Accordingly, the temperature 
and density in the outer region become higher owing to the larger 
plasma heating (Figure 5). 
4.1.2. Wave Polarization 

We examine the effect of the wave polarization. Linearly po- 
larized Alfven waves dissipate by the direct steepening [Hollweg , 
1982; Suzuki , 2004] to MHD fast shocks as well as by the paramet- 
ric decay instability [Goldstein , 1978]. On the other hand, circu- 
larly polarized components do not steepen, and it dissipates only by 
the parametric decay instability in our one-fluid MHD framework, 
so that the dissipation becomes less efficient. 

Figures 7 and 8 compare the results of sinusoidal linearly polar- 
ized (Run IV) and circularly polarized (Run V) fluctuations with 
the same amplitude ((dv±, 0) = 0.7km s _1 ). The results of Run II 
are also plotted for comparison. It is shown that the coronal heating 
and solar wind acceleration are still achieved by the dissipation of 
the circularly polarized Alfven waves (Figure 7) although a large 
fraction of the wave energy remains at the outer boundary (Figure 
8) owing to the less dissipative character as expected. 

A key ingredient for the efficient dissipation is the rapid decrease 
of the density owing to the gravity as discussed in §3. The circularly 
polarized Alfven waves which are initially sinusoidal are quickly 
deformed by the rapid variation of the density in the chromosphere 
and the low corona; simple sinusoidal waves cannot propagate any 
further unless the wavelengths are sufficiently short. The shape de- 
formation indicates that the waves are partially reflected and the 
incoming Alfven waves are generated. Compressive slow waves 
are also easily generated even without the steepening because the 
variation of the magnetic pressure is present owing to the large dif- 
ference of the Alfven speed within one wavelength [Boynton & 
Torkelsson , 1996; Ofman & Davila , 1997, 1998; Grappin et al., 
2002], 

It has been generally believed that the decay instability of the 
Alfven waves is not important in the context of the coronal heating 
and the solar wind acceleration. However, our simulation shows 
that this is not the case; the three-wave interactions are greatly en- 
hanced by the long wavelength effect so that the sufficient wave 
dissipation occurs in the real solar corona. 

The case of the circularly polarized waves gives smaller coronal 
density and temperature since the wave dissipation is less effective. 
However, the solar wind speed is rather faster, up to 700km s _1 at 
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0.1 AU, because the less dense plasma can be efficiently accelerated 
by transfer from the momentum flux of the Alfven waves. 
4.1.3. Wave Amplitude 

We study dependences on the amplitudes of the input fluc- 
tuations at the photosphere. We compare the results of larger 
(dfx.o) = 1.4km s (Run VI) and smaller (cfox.o) = 0.4km 
s _1 (Run VII) cases with the fiducial case ((cfox.o) = 0.7km 
s _1 ; Run II) in Figures 9 and 10. We also show the results of 
(dv ±fi ) = 0.3km s _1 (Run IX) at t = 18.3hr in Figure 9, although 
quasi steady-state behavior is not achieved. At later time in Run IX, 
the density and temperature decrease with time as explained later. 

The maximum temperature of Run VII ((d«x,o) — 0.4km s _1 ) 
is ~ 5 x 10° K, which is cooler than the usual corona. The den- 
sity is much lower than the fiducial case by 1-2 orders of mag- 
nitude because the sufficient mass cannot supply into the corona 
by the chromospheric evaporation owing to the low temperature; 
the evaporation is drastically suppressed as T decreases since the 
conductive flux (F c oc T 5 / 2 ^) sensitively depends on T. As a 
result, the mass flux (pv r ) becomes more than an order of magni- 
tude lower than that of the present solar wind. These tendencies 
are more extreme in Run IX; The coronal temperature and density 
become further low in Run IX; T ~ 2 x 10 5 K and the density in 
the corona and solar wind is smaller by 3 orders of magnitude than 
the fiducial case. 

This behavior can be understood by the wave dissipation (Figure 
10). The input wave energy (oc (dv ± )) in Run VII is 1/3 of that 
of Run II. However, the upper panel of Figure 10 shows that the 
remained S c 's in both cases are similar at 0.1AU, which indicates 
that the Alfven waves do not effectively dissipate in the smaller 
(dv±,o) case. This is because the Alfven speed (= B/ v / 47rp) is 
larger due to the lower density and the nonlinearity, {5va,+)/va, 
becomes weaker. As a result, the final energy flux of the solar wind 

at 0. 1 AU, mostly consisting of the kinetic energy (pv r ^f), becomes 
one order of magnitude smaller although the ratio of the input wave 
energy is only 1/3. 

Once the coronal density starts to decrease, a positive feedback 
operates. Namely, the decrease of the density leads to weaker 
non-linearity of the Alfven waves, which reduces the plasma heat- 
ing by the wave dissipation. This decreases the coronal temper- 
ature, and further reduces the coronal density by the suppression 
of the chromospheric evaporation. This takes place in Run IX 
((dv± : o) = 0.3km s _1 ) so that the coronal density and tempera- 
ture continue to decrease at later time, instead of maintaining steady 
corona and solar wind. 

Our results show that (d«x,o) ~ 0.4km s _1 is the criterion 
of the photospheric fluctuations for the formation of the stable hot 
plasma. To get the maximum coronal temperature ~ 10 6 K, 

{dv±fi} ~ 0.7km s _1 is required. Otherwise if (cfox.o) ~ 
0.3km s~ , the low-frequency Alfven waves cannot maintain the 
hot corona, and the solar wind mass flux becomes drastically small; 
solar winds virtually disappear (see §4.1.4). 

The solar wind speeds are slightly faster in the cases of smaller 
(dv± }0 ) (accordingly smaller energy input), because the densities 
at the coronal bases are much smaller owing to the suppression 
of the chromospheric evaporation; the smaller amounts of the ma- 
terials are accelerated to the faster speeds. The differences be- 
tween the photospheric and coronal base densities are larger in the 
smaller (dvx,o) cases and the amplification of dv± (oc p -1 ^ 4 if 
Alfven waves are nondissipative in static media; e.g. Lamers & 
Cassinelli 1999) is larger. As a result, (dv±) at the coronal base 
is (weakly) anti-correlated with (dux.o) at the photosphere. Thus, 
our result is consistent with the previous calculations [Sandbaek & 
Leer , 1995; Ofman & Davila , 1998] that show the positively cor- 
relation between the Alfven wave amplitude at the coronal base and 
the solar wind speed. 

Larger (d«x,o) gives larger coronal density as shown in Run 
VI. The initial increase of the temperature starts from a deeper lo- 
cation around r ~ 0.005-Rq than the other cases. Thanks to this, 
a decrease of the density is slower (larger pressure scale height) 
so that the density around r = 1.01i?Q is two orders of magni- 
tude larger than that of Run II. However, the temperature decreases 



slightly instead of a monotonical increase; it cannot go over the 
peak of the radiative cooling function at T ~ 10 5 K [Landini & 
Monsignori-Fossi, 1990] because the radiative loss is efficient ow- 
ing to the large density. The second increase of the temperature 
begins from r ~ 1.03-Rq and above there the corona is formed. 
The coronal density and temperature are larger than those in Run 
II. In particular the density in the outer region is 10 times larger, 
and the mass flux of the solar wind is larger by the same extent. 
The top panel of Figure 10 shows that the Alfven waves dissipate 
more effectively in the larger (dfx.o) case by the same mechanism 
as explained above for the smaller (cfox.o) cases. 

Figure 11 summarizes the maximum coronal temperature, 
T max (K) (top), the solar wind mass flux at 0.1AU, (piv)o.iAu(g 
cm _2 s _1 ) (middle), and solar wind speed at 0.1AU, « r ,o.iAu(km 
s _1 ), (bottom) as functions of (cfox,o)(km s _1 ). On the 
right axis of the middle panel, we show the proton flux at 
1AU, (n p u r )iAu(cm _ s _1 ). When deriving (n p iv)iAu from 
(pw r )o.iAu, we use the relation of pv r oc r~ 2 , and assume the 
solar elemental abundance. T max and (pw r )o.iAu have positive 
correlations with (d«x,o)- The dependence of (ptv)o.iAu i s quite 
steep because of the nonlinearity of the Alfven waves as explained 
above. The dependence of T max is more gradual due to the redistri- 
bution of the temperature by the thermal conduction, iv.o.iau has 
weakly negative correlation with (di>x,o) for (cfox.o) > 0.4(km 
s _1 ) because in smaller (d«x,o) cases the coronal density is fur- 
ther lower and smaller amounts of materials are accelerated more 
effectively. For (dv±) = 0.3(km s _1 ), the coronal temperature 
and the mass flux of the solar wind becomes smaller later time and 
the wind speed varies quite a lot; we use upper limits for T max and 
(pVr)o.iAu, and a dashed error bar for iv.o.iau- 

In our simulation we only consider the steady fluctuations with 
constant (cfox,o)'s to study the basic relation between the footpoint 
fluctuations and the properties of the solar winds. In the real situa- 
tion, however, (cfox.o) would vary in time. Our results imply that a 
small temporal variation of (dwx,o) at the surface possibly leads to 
large time dependent behavior of the solar wind. As one example 
we would like to discuss an event of solar wind disappearance in 
the following section. 

4.1.4. Disappearance of Solar Wind 

On May 10-12, 1999, observed solar wind density near the earth 
drastically decreases well below 1 cm -3 for ~ a day and to 0.1 
cm -3 at the lowest level, in comparison with the typical value ~ 5 
cm~ 3 [Le et al., 2000; Smith et al., 2001]. Various mechanisms are 
widely discussed to explain this disappearance event [Usmanov et 
al. , 2000; Richardson et al., 2000; Crooker et al., 2000]. In §4.1.3, 
we have shown that the solar wind density sensitively depends on 
the amplitude of the photospheric fluctuation provided that the non- 
linear dissipation of low-frequency Alfven waves operates in the 
solar wind heating and acceleration. According to Figure 1 1 the 
solar wind density and mass flux becomes ~ 100 and ~ 50 times 
smaller respectively if (dv±,o) changes from 0.7km s _1 to 0.4km 
s _1 ; a small variance of the energy injection at the solar surface 
leads to a large variation of the solar wind density because of the 
nonlinearity. Therefore, we can infer that if the amplitude of the 
photospheric turbulence becomes ~ 1/2 during ~a day, this event 
of the solar wind disappearance is possibly realized. 

The observation shows that the velocity of the sparse plasma dis- 
cussed above is lower than that of the surrounding plasma, which 
is inconsistent with our results (§4.1.3). We are speculating that 
the stream interaction would resolve the inconsistency. The sparse 
region is easily blocked by the preceding dense plasma in the spiral 
magnetic fields (Parker spiral) even if its speed in the inner region 
(say ~ 0.1 AU) is faster [Usmanov et al. , 2000], because the ram 
pressure is much smaller than that of the surrounding plasma. For 
quantitative arguments, however, two dimensional modeling is re- 
quired. 

4.2. Flux Tube properties 



X- 6 



T. K. SUZUKI & S. INUTSUKA: SOLAR WINDS BY ALFVEN WAVES FROM PHOTOSPHERE 



It is widely believed that field strength and geometry of open 
flux tubes are important parameters that control the solar wind 
speed. Wang & Sheeley [1990, 1991] showed that the solar wind 
speed at ~ 1AU is anti-correlated with / max from their long-term 
observations as well as by a simple theoretical model. Ofman & 
Davila [1998] showed this tendency by time-dependent simula- 
tions as well. Fisk et al. [1999] claimed that the wind speed should 
be positively correlated with B r x) by a simple energetics consider- 
ation. Kojima et al. [2005] have found that the solar wind velocity 
is better correlated with the combination of these two parameters, 
B r .o//max, than l// max or B r fi from the comparison of the out- 
flow speed obtained by their interplanetary scintillation measure- 
ments with observed photospheric field strength. Suzuki [2004] 
and Suzuki [2006] also pointed out that B rj0 / /max should be the 
best control parameter provided that the Alfven waves play a dom- 
inate role in the coronal heating and the solar wind acceleration. 

Here we test the last scenario by our simulation. In Figures 
12 and 13, we compare the results of a flux tube with smaller 
Bro//max = 322(G)/450 (Run VIII) with the fiducial case 
(Br,o//max = 161(G) /75) (Run II) for the polar coronal hole. We 
chose the parameters of Run VIII with equatorial and mid-latitude 
coronal holes in mind; they are surrounded by closed structures, 
and thus, have moderately larger photospheric field strength and 
much larger flux tube divergence than the polar coronal hole. 

Figure 12 shows that Run VIII gives much slower solar wind 
speed which is consistent with Kojima et al. [2005], Suzuki [2004], 
and Suzuki [2006], with slightly hotter and denser corona. These 
results can be understood by positions of the wave dissipation. The 
upper panel of Figure 13 indicates that the outgoing Alfven waves 
dissipate more rapidly in the smaller B r ,o//max case. This is be- 
cause the nonlinearity of the Alfven waves, (Sva,+)/va is larger 
due to smaller v\ oc B r (oc B r ,o// ma x in the outer region where 
the flux tube is already super-radially open) even though the abso- 
lute amplitude {{dv±}) is smaller (bottom panel of Figure 12). As 
a result, more wave energy dissipate in the smaller B r s> / /max case 
in the subsonic region and less energy remains in the supersonic 
region. In general, energy and momentum inputs in the supersonic 
region gives higher wind speed, while those in the subsonic region 
raises the mass flux (pv r ) of the wind by an increase of the density 
[Lamers & Cassinelli, 1999]. Therefore, the smaller B r ,o// m ax 
case gives slower wind with higher coronal density, whereas the 
solar wind density in the outer region is similar to that of the larger 
B r ,o//max case on account of the dilution of the plasma in the 
more rapidly expanding flux tube. We can conclude that B r fi / /max 
controls the solar wind speed and the coronal density through the 
nonlinear dissipation of the Alfven waves. 

4.3. Summary of Parameter Studies 

We have examined the dependences of the coronal and wind 
properties on the various wave and flux tube parameters in §4.1 
and §4.2. One of the important results is that we do not find any 
subsonic 'breeze' solution in our simulations. All the Runs with 
(dv±fi) ~ 0.3km s _1 show the transonic feature, hence, tran- 
sonic solar winds are natural consequences of the dissipation of 
the low-frequency Alfven waves. Smaller energy injection only re- 
duces the density (and consequently the mass flux) of the wind. 
The wind speed itself is rather slightly faster for smaller wave en- 
ergy because the coronal base density becomes lower owing to the 
suppression of the chromospheric evaporation. 

From §4.1.1 and §4.1.2, we can conclude that the dependences 
on the spectra and polarizations of the input fluctuations are weak 

as long as the low-frequency ( ~ 0.05Hz) waves are considered. 
(dv±fi) and -B ri o//max are the important control parameters, and 
the dependences are summarized below: 

• The coronal temperature and the density in the corona and 
solar wind are mainly determined by the photospheric amplitude, 

(dv±,o). The corona with temperature ~ 10 6 K is formed if 

(dt>x,o) ~ 0.7km s _1 . Larger (cfox.o) gives higher density in 

the corona and the solar wind. If (cfox.o) ~ 0.3km s _1 the hot 
plasma cannot be maintained and the mass flux of the solar wind is 
unrealistically small. 

• The solar wind speed is mainly controlled by Br,o//ma X ; 
faster winds come from open flux tubes with larger B r fi //max- 



5. Fast and Slow Solar Winds 

The observed solar winds can be categorized into two distinctive 
types. One is the fast solar wind which is mainly from polar coro- 
nal holes and the other is the slow wind which is from mid- to low- 
latitude regions. Apart from the velocity difference, both density 
of solar winds [Phillips et al., 1995] and freezing-in temperatures 
of heavy ions [Geiss et al., 1995] which reflect the coronal tem- 
peratures are anti-correlated with velocities. Origins of the slow 
solar wind are still in debate; there are mainly two types argued 
(see Wang et al.2000 for review). One is the acceleration due to 
intermittent break-ups of the cusp-shaped closed fields in the equa- 
torial region (e.g. Endeve, Leer, & Holtzer 2003). The other is the 
acceleration in open flux tubes with large areal expansions in low- 
and mid-latitude regions. Kojima et al. [1999] detected low-speed 
winds with single magnetic polarities, originating from open struc- 
ture regions located near active regions, which indicates a certain 
fraction of the slow stream is coming from open regions, similarly 
to the fast wind. In this paper we focus on the latter case with 
respect to the slow solar wind. 

Based on our parameter studies summarized in §4.3, we can in- 
fer the suitable (cfox.o) and i? r ,o//max for the slow solar wind by 
comparing with those for the fast wind (Run 1/11; also in paper I); 
the slow wind requires larger (di>x,o) and smaller B r ,o //max than 
the fast wind. The adopted parameters are summarized in table 1 
(Run X). 

In the panels on the left side of Figure 14 we compare the results 
of Run X (slow wind) and Run I (fast wind) overlayed with recent 
observations of slow solar winds (see Figure 2 for the observations 
of fast winds). The temperature and density of the slow wind case 
becomes larger on account of the larger (cfox,o), which explains the 
observations. On the other hand, smaller Bo// m „ gives slower 
terminal speed. As a result, the observed anti-correlation of the 
wind speed and the coronal temperature [Schwadron & McComas 
, 2003] is well-explained by our simulations. In the slow wind case 
(Run X), the acceleration of the outflow is more gradual, and it is 
not negligible in r ~ 20© (e.g. Nakagawa et al.2005). 

In the panels on the right of Figure 14, we show the properties 
of the wave dissipation. The top right panel compares transverse 
amplitude, (dv±), averaged over 28 minutes in the fast and slow 
wind cases. (dv±) is larger in the fast wind. This is because the 
Alfven waves become less dissipative in the fast wind conditions; 
larger B(oc Bo//max) and smaller p give larger va so that the 
nonlinearity of the Alfven waves, {8va,+)/va, is systematically 
smaller. As a result, the Alfven wave dissipates more slowly in 
the fast wind case, which is shown in the plot of S c of the wave 
(the bottom right panel of Figure 14). This tendency of the larger 
(dv±) in the fast wind is expected to continue to the outer region, 
~ 1 AU. Various in situ observations in the solar wind plasma near 
the earth also show that the fast wind contains more Alfvenic wave 
components (see Tsurutani & Ho 1999 for review), which can be 
explained via the nonlinear dissipation of Alfven waves based on 
our simulations. 

The middle right panel of Figure 14 shows longitudinal fluc- 
tuation, {dvu), averaged over 28 minutes, (dvu) reflects the am- 
plitude of slow MHD waves generated from the Alfven waves via 
the nonlinear process (§3). Slow waves are identified in polar re- 
gions [Ofman et al., 1999] and low-latitude regions [Sakurai et al., 
2002] at low altitudes (r ~ 1.2J?q). The observed amplitudes 
are still small (~ 7.5km s _1 in the polar regions and ~ 0.3km 
~ 1 in the low-latitude regions) because of the large density there, 
which is consistent with our results. What is more interesting is 
longitudinal fluctuations in outer regions, our results exhibits that 
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{dv\\ } in the solar wind plasma is not small. Typically, the simu- 
lations give {dv\\ } ~ lOOkm/s in 37?© ~ r ~ 107?© of the 

fast wind and {dv\\) ~ 20km/s in 2R Q ~ r ~ lOR© of the 
slow wind, whereas these might be modified if we take into account 
multi-dimensional effects (§6.2). This is directly testable by in situ 
measurements of future missions, Solar Orbiter and Solar Probe, 
which will approach to ~45 and 4 R Q , respectively, corresponding 
to the inside of our computation domain. 

Our simulations show that the different types of the solar winds 
can be explained by the single process, the dissipation of the low- 
frequency Alfven waves, although we do not intend to exclude 
other possibilities. The varieties of the solar winds are due to the 
varieties of the footpoint amplitudes as well as the magnetic fields 
and geometry of the flux tubes. Our choice of -B ri o//max is con- 
sistent with the obtained data by Kojima et al. [2005] who report 
that the slow winds are mainly from the low-to mid-latitude coronal 
holes with smaller B r ,o//max. In contrast, quantitative arguments 
on (dv± } o) need more detailed fine scale observations of the mag- 
netic fields with high cadence by future telescopes such as Solar-B. 

6. Discussions 

6.1. Wave Generation 

In this paper we have assumed that the origin of the waves is 
the steady turbulent motions at the photosphere. However, tran- 
sient activities are also expected to play a role in the wave genera- 
tion. Sturrock [1999] proposed small flare-like events triggered by 
magnetic reconnections of closed loops in the chromosphere ex- 
cite MHD waves in the corona. Miyagoshi et al. [2005] also show 
that interactions between open field lines and emergent flux tubes 
excite Alfven waves at a location above the photosphere. These 
waves could directly heat up the corona since they do not suffer the 
attenuation in the chromosphere and the TR. 

An important issue for the wave generation by these transient 
events is the energetics. The energy release from each event is gen- 
erally thought to be small, categorized as a micro- or nano-fiare 
[Parnell & Jupp, 2000; Katsukawa & Tsuneta, 2001]. To clarify 
how these waves dominantly work in the heating and accelera- 
tion of the solar wind, we should determine the frequency of these 
events. These small-scale events might be important in the heating 
of not only the open coronal holes but the closed regions. The de- 
termination of the total energy release from small transient events 
by future observations is quite important in terms of the coronal 
heating in various portions of the corona. 

6.2. Limitation of Our Simulations 

We have shown by the self-consistent simulations that the dissi- 
pation of the low-frequency Alfven waves through the generation 
of the compressive waves and shocks is one of the solutions for 
the heating and acceleration of the plasma in the coronal holes. 
However, the validity of the ID MHD approximation needs to be 
examined. 

We think that the MHD approximation is appropriate as a whole 
for studies of the average properties of the plasma. Let us esti- 
mate the Larmor radius, l Lalm or, of the protons which compose the 
main part of the plasma due to fluctuations of field lines, B±, un- 
der the coronal condition. Typically, B± ~ 0.lB r ~ 0.1G, and 
then ?Larmor ~ 0.1km for a thermal proton with ~ 100km s _1 . On 
the other hand the shortest wavelength, A m i n , of the Alfven waves 
we are considering is A m i n ~ 3 x 10 4 km (a period of 20s and 
va ~ 1500km s _1 .) and the simulation grid size in the corona 
is ~ 3 x 10 km (~ 1/10 of A m i n ); the thermal protons are sup- 
posed to be well randomized through turbulent magnetic fields on 
the scales we are dealing with. 

In our simulations the heating of the plasma is done by MHD 
shocks (§2). Our simulations cannot resolve detailed structures 
of shocks since the size of the grid is larger than the width of 
the shock front which is an order of a Larmor radius. However, 



the global structure of shocks can be appropriately treated by our 
simulations. The thermal particles could be 'trapped' around the 
shock regions by the random magnetic fields so that the global 
properties of shocks would satisfy the MHD condition. Therefore, 
the shock heating rate calculated in our MHD simulations is sup- 
posed to give a reasonable estimate in the coronal region. In the 
outer region (r ~ 2Rq), /Larmor (oc 1/B) increases and A m i n 
(~ t>A/fmax oc oc B / \fp) decreases as r increases. Even 
around 0.3AU, tarmor is about ~1/10 of A m i n ; the approximation 
still remains reasonable. 

However, we should be cautious about the fact that the shocks 
are collisionless; mean free path due to Coulomb collisions, which 
is ~ 100km for electron-electron collisions and even larger for 
electron-proton and proton-proton collisions under the typical coro- 
nal condition (T = 10 6 K and n = 10 9 cm ), is much larger 
than ^Larmor so that the particles are not thermalized by Coulomb 
collisions but only randomized by fluctuating fields. In such a 
condition, the particle distribution function possibly deviates from 
the Maxwell-Boltzmann distribution, which is actually observed 
in solar wind plasma [Marsch et al., 1982a]. Non-thermal, or 
even supra-thermal, particles might modify the heating rate, if their 
amount is not negligible. The shock energy is transfered not only 
to the thermal heating but also to acceleration of the nonthermal 
component. Although this is beyond the scope of the present paper, 
this issue should be carefully considered in more detailed models. 

Effects of multi-component plasma also need to be taken into 
account [Lie-Svendsen , 2001; Of man , 2004]. The thermal con- 
duction and radiative cooling are mainly involved with electrons, 
and the electron temperature becomes systematically lower than 
the ion temperature without effective thermal coupling by Coulomb 
collisions. The observed proton temperature is actually higher than 
the electron temperature even around a few solar radii [Esser et al., 
1999]. The temperature in our simulations represents the electron 
temperature because it is mainly determined by the thermal con- 
duction and radiative cooling. The actual thermal pressure is sup- 
posed to be higher than the simulated temperature due to the higher 
proton temperature. Therefore, the acceleration profile might be 
slightly modified whereas the the dynamics are mostly controlled 
by the magnetic pressure of the Alfven waves. 

Collisionless processes would also be important for the wave 
dissipation [Tsurutani et al., 2005]. The heating rates by the 
dissipation of Alfven waves are different for electrons and ions 
[Hasegawa & Lui , 1975; Tsiklauri et al., 2005]. The compres- 
sive waves generated from the Alfven waves would suffer transit- 
time damping by interactions between magnetic mirrors and surfing 
particles [Barnes, 1966; Suzuki , 2006]. If high-frequency (ioncy- 
clotron) waves are produced by frequency cascade, the heating sen- 
sitively depends on mass-to-charge ratios of particles [Dusenbery 
& Hollweg , 1981; Marsch et al., 1982b]. It is believed that this 
process can explain the observed preferential heating of the per- 
pendicular temperature of the heavy ions [Kohl et al., 1998], while 
there are potential difficulties in the proton heating as discussed in 
§1- 

In the chromosphere, the gas is not fully ionized. In such cir- 
cumstances, the friction due to ion-neutral collisions might be im- 
portant in the dissipation of Alfven waves [Lazarian et al., 2004; 
Yan & Lazarian , 2004]. 

Generally, the shock dissipation tends to be overestimated in ID 
simulations because the waves cannot be diluted by the geometri- 
cal expansion. On the other hand, there are other mechanisms of 
the wave damping due to the multi-dimensionality [Ofman , 2004], 
such as phase mixing [Heyvaerts & Priest, 1983; Nakariakov et 
al., 1998; De Moortel et al., 2000] and refraction [Bogdan et al., 
2003]. Moreover, the solar wind plasma is more or less turbulent 
and it must be important to take into account the plasma heating 
through cascades of Alfvenic turbulences in the transverse direction 
[Goldreich & Sridhar , 1995; Mattaeus et al., 1999; Oughton et al., 
2001]. Studies of the imbalanced cascade, in which a wave com- 
ponent in one direction has larger energy than that in the counter 
direction, seem to be important for us to understand transport phe- 
nomena in the solar wind plasma [Lithwick & Goldreich , 2003], 
though it has not been fully understood yet. 
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Therefore, the Alfven waves might also be dissipated by mech- 
anisms different from those included in our simulations. Accord- 
ingly, variation of wave amplitudes might be modified when these 
additional dissipation processes are considered. Self-consistent 
simulations including these various processes remain to be done 
in order to arrive at final conclusion. 

7. Summary 

We have performed parametric studies of the coronal heating 
and solar wind acceleration by the low-frequency Alfven waves in 
open coronal holes. We have performed ID MHD numerical sim- 
ulations from the photosphere to 0.3 or 0.1 AU. The low-frequency 
Alfven waves with various spectra, polarization, and amplitude are 
generated by the footpoint fluctuations of the magnetic field lines. 
We have treated the wave propagation and dissipation, and the heat- 
ing and acceleration of the plasma in a self-consistent manner. 

We have found that the transonic solar winds are accomplished 
in all the simulation Runs. Smaller energy injection does not re- 
duce the outflow speed but the density to maintain the transonic 
feature. The atmosphere is heated up to the coronal temperature 

( ~ 10 6 K) if the photospheric amplitude, (dvx.o) ~ 0.7km s _1 . 
Otherwise, the temperature and density become much lower than 
the present coronal values. If (dv± t o) ~ 0.3km s _1 , the temper- 
ature becomes less than a few 10 5 K and the sufficient mass cannot 
be supplied into the corona owing to the suppression of the chro- 
mospheric evaporation. The stable hot corona cannot be maintained 
any longer, and the mass flux of the solar wind becomes at least 3 
orders of magnitude smaller than the observed value of the present 
solar wind. This shows that the solar wind almost disappears only 
by reducing the photospheric fluctuation amplitude by half. On 
the other hand, the case with (d«x,o) = 1.4km s~ gives 10 times 
larger density than the fiducial case ((dt>x,o) = 0.7km s _1 ). These 
sensitive behaviors of the solar wind mass flux on (dvx.o) can be 
explained via the nonlinear dissipation of Alfven waves. Our simu- 
lations have also confirmed that the positive correlation of the solar 
wind speed with B r ,o//max obtained by Kojima et al. [2005]. 

We have finally pointed out that both fast and slow solar 
winds can be explained even solely by the dissipation of the low- 
frequency Alfven waves with different (dt>x,o) and B r ,o//max, 
whereas we do not intend to exclude other possibilities. Fast winds 
are from flux tubes with larger B r ,o/fmax and smaller (dv±,o), 
while slow winds are from flux tubes with smaller B r ,o/fmax 
and larger (dv±,o). These choices naturally explain the follow- 
ing observed tendencies: (i) the anti-correlation of the solar wind 
speed and the coronal temperature, and (ii) the larger amplitude of 
Alfven waves in the fast wind in the interplanetary space. The ten- 
dency with respect to B r ,o//max is consistent with the observed 
trend [Kojima et al., 2005]. To determine (cfox.o) in various coro- 
nal holes, fine scale observations of magnetic fields with high ca- 
dence are required; this is one of the suitable targets for Solar-B 
which is to be launched in 2006. 
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where s = p/p 1 and c s = y/jp/p, and subscripts, _L, 1 and _L, 2, 
denote first and second transverse components. C consists of terms 
due to external force, additional heating and cooling, and curvature. 
In our case, we have taken into account gravity, thermal conduction, 
radiative cooling, and superradial expansion of flux tubes: 



B r B I o \ 



Fl(v r V±,l 

Fi{v r v ± , 2 j^r- ) 

Fi(B± A v r - B r «i,i) 
Fi(B±, 2 v r - B r v±, 2 ) 

F 2 pV r 



where Fi = -J-=-§-( r ^f) and F 2 = i£(r 2 /) are aris 



ry/f 

ing from the geometrical effect, and B\ 

2 2 i 2 

V± = «L,1 + «x,2 

the thermal conduction term (— 

_pr A f or 



B , 



,1 + B ±,2 



and 



Please note that C is negligible except for 
(r 2 /F c )) in our simulations 
if r ou t is set at a sufficiently distant location. 

Using eigen values, \i (i = 1, • • • , 7), and eigen vectors, Ii, 
defined as 



IiA = Xili, 



(A2) 



we can derive characteristic equations for the seven MHD waves 
from eq.(Al) : 

I i ^-+X i li^+I i -C = (A3) 
at or 

Xi corresponds to phase speed of each wave, 

Al = V r + Vf, A 2 = V r + V A , A 3 = V r + V s , A 4 = V r , 

A 5 = V r — V s , X 6 = V r — VA, X 7 = V r — Vf, (A4) 
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We introduce our prescription for the outgoing boundary con- 
dition of the MHD waves. The method is similar to that in Wu et 
al. [2001], while ours is implemented for stable treatment of strong 
MHD shocks. 

Basic MHD equations can be expressed in a matrix form as 



7j=> and v s = ^l^ + ^-p- V^+i^) 2 -^] are 
phase speeds of fast, Alfven , and slow modes, respectively. Here 
we define total magnetic field by B 2 = B 2 + B\ . The eigen vec- 
tors are explicitly written as 



0. 
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The characteristic equations enable us to implement the outgo- 
ing boundary condition in a clear and direct manner. All the incom- 
ing waves should be removed for the proper outgoing boundary. 
This is accomplished by setting physical variables, U, to be con- 
stant in space for the incoming characteristics [Thompson , 1987]. 
Then, we simply impose following condition at r = r out : 



0, 



where 



U = 



\ . r. OU 



A t > (Outgoing) 
A, < (Incoming) 



(A5) 



(A6) 



In our scheme, we firstly perform time evolution to derive 
U™(1 < j < Jout — 1) at time n except the outermost grid, j out , 
from E/" _1 (l < j < jout) at the previous time step, n — 1, by the 
second order MHD Godunov-MOCCT scheme. The physical vari- 
ables at the outermost grid, U" out are determined at the end of each 
time step by using C7™~ and U" out _ 1 - At this time we integrate 
eq.(A5) implicitly in time for numerical stability: 



At 



0. 



(A7) 



We adopt upwind discretization for the spatial derivative ap- 
peared in d for the outgoing characteristic, (^r)" out = 

t ^"" r A ?"" t ' 1 . Note that (7i)" ut - (£i)£ at - and C™ ut are func- 
tions of Uj" ont , so we need iteration to determine U" out ■ Equation 
(A7) determines the physical quantities at the outermost mesh to 
satisfy the outgoing condition for all the seven MHD characteris- 
tics. 

We carry out reflection tests for our implementation by using 
different types of waves with various amplitudes from linear to ex- 
tremely nonlinear regimes. We here consider two types of waves 
as typical examples : (i)Alfven waves, and (ii) oblique fast waves 
which propagate in 45 degree with respect to the underlying mag- 
netic field. We use Cartesian coordinates with homogeneous back- 
ground for the tests. Solitary waves which travel in the right direc- 
tion are initially set, and the outgoing condition is implemented at 
the right boundary. In order to give the waves which purely prop- 
agate in one direction even in the extremely nonlinear regime, we 
adopt simple wave solutions. Circularly polarized Alfven waves 



traveling in one direction are the exact solutions of the basic MHD 
equations. For the oblique fast waves, p, p, v x , and v y can be ex- 
pressed by By [Wu, 1987]: 
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(All) 



Note that we can treat fast (and slow) waves in two dimensional 
(x&zy) space without loss of generality because they are coplanar, 
while the three dimensional space is required for circularly polar- 
ized Alfven waves. We initially give sinusoidal variations of B y for 
the fast waves and the other quantities are determined by Equations 
(A8) - (All). If the wave amplitude, 5v, is sufficiently smaller than 
Vf, all the quantities can be approximated by sinusoidal variations. 

Figure 15 shows an example of the reflection tests for the 
oblique fast waves. We non-dimensionalize the B values by 
1/4-7T — > 1. The initial background conditions are, p = 1, p = 1, 
B x — 1.5, By = 1.5, v x — 0, and v y — 0, which give plasma 
/3 value, /3 = 2p/B 2 (= 8np/B 2 in cgs-Gauss unit)= 0.44. We 
give magnetic field amplitude, SB y = 6.5, which corresponds to 
relatively strong non-linearity, Sv/vf = 3.1. The simulation box is 
from x = to x = 4 and the number of grid points is 512. The 
initial solitary wave is put between x — 2 and 3. The initial condi- 
tion (t — 0) and results at t = 0.6 are plotted in dashed and solid 
lines. To estimate errors at the right boundary, we also carry out 
simulation in larger box between x = and 8(dotted lines), which 
is free from the boundary errors at x = 4. 

Because of the strong nonlinearity, the wave rapidly steepens to 
form shocks before reaching x = 4. These shocks travel to the 
right direction, which are seen in the simulation with the larger box 
size (dotted lines). At t = 0.6, the main part of the wave goes 
out of the simulation box (x = 4) except the trailing edge. The 
differences between solid and dotted lines indicate errors due to 
our implementation of the outflow boundary condition. Figure 15 
shows the errors are quite small; relative errors to the initial ampli- 
tude are at most 1% for each variable. Our implementation is good 
enough for the passage of the moderately strong shocks. 

Figure 16 shows the relative errors to the initial amplitude, 
5v alIOI /Sv on the normalized wave amplitude, Sv/vf and Sv/v A , 
for the oblique fast waves and circularly polarized Alfven waves. 
As for the relative errors, we plot the errors of v y , whereas the 

other quantities (p,p, B y ,v x ) give similar results. If Sv/vp-^ ~ 3 
(v p h = Vf or v A is phase speed), the relative errors are smaller 
than 0.01 and our implementation is sufficiently acceptable. The 
Alfven wave with 8v/v A = 10 gives larger 8v CIIOI /Sv ~ 0.05. 
This is because as the wave travels it dose not only consist of 
Alfven modes but contains non-linear fast and slow wave compo- 
nents. However, this value (8v CIIOI /8v ~ 0.05) is still acceptable 
since an error in energy (oc Sv 2 ) is only 0.3%. 

Notes 

1 . It is instructive to note that the incoming Alfven waves have the positive 
correlation between v± and B ± (dark-dark or light-light in the figures), 
while the outgoing modes have the negative correlation (dark-light or 
light-dark). 

2. The phase correlation of the longitudinal slow waves is opposite to that 
of the transverse Alfven waves. The outgoing slow modes have the pos- 
itive correlation between amplitudes of v r and p, (8v r Sp > 0), while 
the incoming modes have the negative correlation (8v r Sp < 0). 
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Figure 1. The left figure shows the geometry of the open flux 
tube in our simulations. Arch represents the solar surface and 
solid lines indicates magnetic field lines. Dashed lines corre- 
spond to radial expansion. The flux tube expands super-radially 
around Ri and R2. This geometry mimics realistic open flux 
tubes on the sun. An example is shown on the right. A siz- 
able fraction of the surface is occupied by closed loops with a 
typical length ~ 10 4 km (=O.O14R ). Above that height open 
structures become more dominant. Then, the open flux tube ex- 
pands super-radially at ~ Ri(= 1.01i?Q in our simulations). 
The tube further expands at ~ R2(~ 1.2Rq in our simulations) 
due to the large-scale dipole structure. 
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Figure 2. Results of Run I with observations of fast solar 
wind. From top to bottom, outflow speed, iv(km s -1 ), tem- 
perature, T(K), density in logarithmic scale, log(p(g cm -3 )), 
and rms transverse amplitude, (dv±)(km s -1 ) are plotted. 
Observational data in the third panel are electron density, 
\og(N e (cm~' i )) which is to be referred to the right axis. 
Dashed lines indicate the initial conditions and solid lines are 
the results at t = 2573 minutes. In the bottom panel, the ini- 
tial value ((dt>±) = 0) dose not appear, first: Green vertical 
error bars are proton outflow speeds in an interplume region by 
UVCS/SoHO [Teriaca et al., 2003]. Dark blue vertical error 
bars are proton outflow speeds by the Doppler dimming tech- 
nique using UVCS/SoHO data [Zangrilli et al., 2002]. A dark 
blue open square with errors is velocity by IPS measurements 
averaged in 0.13 - 0.3AU of high-latitude regions [Kojima et 
al. , 2004]. Light blue data are taken from Grail et al. [1996]; 
crossed bars are IPS measurements by EISCAT, crossed bars 
with open circles are by VLBA measurements, and vertical error 
bars with open circles are data based on observation by SPAR- 
TAN 201-01 [Habbaletal., 1994], second: Pink circles are elec- 
tron temperatures by CDS/SoHO [Fludra et al., 1999]. third: 
Circles and stars are observations by SUMER/SoHO [Wilhelm 
et al., 1998] and by CDS/SoHO [Teriaca et al., 2003], respec- 
tively. Triangles [Teriaca et al., 2003] and squares [Lamy et 
al., 1997] are observations by LASCO/SoHO./oarf/z: Blue cir- 
cles are non-thermal broadening inferred from SUMER/SoHO 
measurements [Banerjee et al., 1998]. Cross hatched region 
is an empirical constraint of non-thermal broadening based on 
UVCS/SoHO observation [Esser et al., 1999]. Green error bars 
are transverse velocity fluctuations derived from IPS measure- 
ments by EISCAT[Canals et al, 2002]. 
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Figure 3. Energy flux (top) and momentum flux (bottom) as a 
function of distance from the photosphere. Each component is 
averaged with respect to time during 28min (the longest wave 
period considered). Each value is normalized by cross section, 
A c , of the flow tube at r = r c (= l.O27? ); note that the 
real flux is larger (smaller) in r < r c (r > r c ). top: Solid, 
dashed, dot-dashed, and dotted lines denote the energy flux of 
the Alfven wave, the enthalpy flux, the kinetic energy flux, and 
the integrated radiative loss, bottom Solid and dashed lines in- 
dicate the advected pressure terms of the Alfven waves, (paiv) 
and the thermal pressure, (pv r ). 
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Figure 4. r — t diagrams for v r (upper-left), p (lower-left), v± 
(upper-right), and B± / B r (lower-right.) The horizontal axises 
cover from Rq to 15Rq, and the vertical axises cover from 
t = 2570 minutes to 2600 minutes. Dark and light shaded 
regions indicate positive and negative amplitudes which ex- 
ceed certain thresholds. The thresholds are dv r = ±96km/s 
for v r , dp/ p = ±0.25 for p, v± = ±180km/s for v±, and 
Bj_/B r = ±0.16 for B±/B r , where dp and dv r are differ- 
ences from the averaged p and v r . Arrows on the top panels in- 
dicate characteristics of Alfven , slow MHD and entropy waves 
at the respective locations (see text). 
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Figure 5. Structures of corona and solar wind for various fluc- 
tuation spectrums at the photosphere. From top to bottom, we 
plot solar wind speed, v r (km s _1 ), temperature, T(K), density 
in logarithmic scale, log(p(g cm' 3 )), and rms transverse ve- 
locity, (dv±)(km s _1 ). Dotted, solid, and dashed lines are re- 
sults of P(y) oc v~ x (Run II), v° (Run III), and 8(1 /v - 180s) 
(Run IV), respectively. Each variable is averaged with respect 
to time during 28min (the longest wave peiriod considered). 
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Figure 6. Adiabatic constant, S'cCerg cm _2 r L ), (top) and 
normalized amplitude, {5va,+)/va, (bottom) of outgoing 
Alfven waves for various spectrums. Dotted, solid, and dashed 
lines are results of P(v) oc v~ x (Run II), v° (Run III), and 
5(l/v — 180s) (Run IV), respectively. Each variable is averaged 
with respect to time during 28min (the longest wave peiriod 
considered). 
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Figure 7. Same as Figure 5 but for different polarizations of the 
perturbations. Solid, dashed, and dotted lines are results of cir- 
cularly polarized sinusoidal waves (Run V), linearly polarized 
sinusoidal waves (Run IV), and linearly polarized waves with 
power spectrum, P(y) oc v^ 1 (Run II), respectively. 
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Figure 8. Same as Figure 6 but for different polarizations of 
the perturbations. Solid, dashed, and dotted lines are results of 
circular polarized sinusoidal waves (Run V), linearly polarized 
sinusoidal waves (Run IV), and linearly polarized waves with 
power spectrum, P{v) oc v^ 1 (Run II), respectively. 
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Figure 10. Same as Figure 6 but for different (dv± } o). Solid, 
dotted, and dashed lines are results of (dv± t o) = 1-4 (Run VI), 
0.7 (Run II), and 0.4(km s _1 ) (Run VII), respectively. 
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Figure 11. The maximum temperature, T max (K), (top), the 
mass flux of solar wind at 0.1AU, (ptv)o.iAu(g cm _2 s _1 ), 
(middle), and the solar wind speed at 0.1AU, uo.iAu(km s _1 ) 
(bottom). On the right axis of the middle panel, we show pro- 
ton flux at 1AU, (n p iv)iAu (cm~ s~ ) estimated by assuming 
the solar elemental abundance. The shaded region in the middle 
panel is observed proton flux at 1 AU. The upper limits in the top 
and middle panels for (dv±fi) = 0.3km s _1 indicates that sta- 
ble corona with the sufficient mass supply cannot be maintained 
and both T max and (piv)o.iAu decreases as the simulation pro- 
ceeds. The wind speed also varies a lot, hence, the dashed error 
bar is used only for that case. 
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Figure 13. Same as Figure 6 but for different B r< o // max . Solid 
and dashed lines are results of -Br-,o(G)// max = 322/450 (Run 
VIII) and 161/75 (Run II), respectively. 
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Figure 14. Simulation results of the slow (Run X; solid lines) 
and fast (Run I; dashed lines) solar wind models in comparison 
with slow-wind observations. On the left from top to bottom, 
outflow velocity, v r (km s _1 ), temperature, T(K), and electron 
density, iV^cm -3 ), are plotted with observation of mid- to low- 
latitude regions where the slow wind comes from, (see Fig- 
ure 2 for the fast-wind observations). On the right from top 
to bottom rms transverse velocity, {dv±)(km s" 1 ), rms lon- 
gitudinal velocity, (dvu)(km s _1 ), and the adiabatic constant, 
Scterg cm~ 2 s~ 1 ) of the outgoing Alfven waves are plotted. 
Observational data; top-left : Shaded region is observational 
data in the streamer belt [Sheeley et al., 1997]. Middle-left : 
Filled diamonds and triangles are electron temperature obtained 
from the line ratio of Fe XIII/X in the mid-latitude streamer 
by CDS/SOHO and UVCS/SOHO respectively [Parenti et al., 
2000] bottom-left : Filled diamonds and triangles are data re- 
spectively from CDS and UVCS on SOHO observation of the 
mid-latitude streamer [Parenti et al., 2000], and filled pentagons 
derived from observation of the total brightness in the equator 
region by LASCO/SOHO [Hayes et al., 2001]. 
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Figure 15. Reflection test for the oblique simple fast wave with 
Sv/vf = 3.1. From top to bottom, v x , v y , B y , p, and p are 
plotted. In the panels on the right each quantity on the left is 
zoomed-up. Dashed and solid lines denote the initial and final 
(t = 0.6) states. Dotted lines are results of simulation in the 
broader region at t = 0.6. Differences between solid and dotted 
lines indicate errors due to our implementation of the outgoing 
boundary condition. 
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Figure 16. Reflection test for the simple fast and the 
Alfven waves. Horizontal axis shows the input wave ampli- 
tudes normalized by the phase speed. Vertical axis plots errors 
due to our implementation of the outgoing boundary condition 
normalized by the input amplitude, (see text) 



